A Lossy Bosonic Quantum Channel with Non-Markovian Memory 



Oleg V. Pilyavets, 1 ' 2 ^ Vadim G. Zborovskii, 3 '0 and Stefano Mancini 1 '^ 

Dipartimento di Fisica, Universita di Camerino, 1-62032 Camerino, Italy 
2 P. N. Lebedev Physical Institute, Leninskn Prospect 53, Moscow 119991, Russia 
3 Troitsk Institute for Innovation and Fusion Research, Troitsk 142190, Russia 

We provide a simple and realistic model to study memory effects in a lossy bosonic quantum 
channel over arbitrary number of uses. The noise correlation among different uses is introduced by 
contiguous modes interactions which results in an exponential decay of the correlations over channel 
uses (modes). This model allows us to characterize the asymptotic behavior of the channel for 
classical information transmission. 
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(\) • I. INTRODUCTION 

£SJ ' A mcmoryless quantum communication channel makes the fundamental assumption that the noise between con- 
secutive uses of the channel is independent. This assumption is reasonable for many real-world applications, but for 
many others the noise may be strongly correlated between uses of the channel. Recently a lot of efforts have been 
dedicated to the development of quantum models that encompass memory effects (see [l[ for an overview). The main 
O ■ motivation that has led to investigate such effects in quantum channels has been the possibility to enhance their 
classical capacity by means of entangled inputs. Such a possibility has been recently put forward in channels with 
continuous alphabet [3, S 0, S @| • These make use of bosonic field modes whose phase space quadratures enable for 
continuous variable encoding/decoding However, most of these works were limited to two channel uses in order to 
provide a proof of principle of the behavior of quantum memory channels. Since the notion of capacity is intimately 
i related with the asymptotic behavior of a channel, there is a persistent wish to move on from small to large (towards 

infinite) number of channel uses. 
t— I ■ The lossy bosonic channel, which consists of a collection of bosonic modes that lose energy en route from the 
transmitter to the receiver, belongs to the class of Gaussian channels which provide a fertile testing ground for the 
general theory of quantum channels' capacities @ and are easy to implement experimentally 0. Here we present a 
' model for a lossy bosonic quantum channel with memory that can be employed for an arbitrary number of channel 
, uses. The memory effects are realized by considering quantum correlations among environments acting on different 
channel uses • This model allows us to characterize the asymptotic behavior of the channel for classical information 
transmission. We then show the usefulness of entangled inputs in presence of memory. In particular, we show that they 
can enhance the classical capacity above that of the corresponding memoryless channel. We also show the utility of 
entangled inputs for rates achievable by conventional decoding procedures (heterodyne and homodyne measurements), 
• • ■ although these do not exceed the classical capacity of the corresponding memoryless channel. 

The paper is organized as follows. In Section II we briefly review the lossy bosonic channels. In Section III we 
introduce the model for the memory effects. Then, Section IV is devoted to the study of the Holevo-x quantity, while 
Section V concerns rates achievable with conventional decoding procedures. Finally, Section VI is for conclusions. 
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II. LOSSY BOSONIC CHANNELS 

The model for a lossy bosonic channel C is depicted in FigJT] On each use of the channel we have an input mode. 
Let us consider n of such modes (uses) with corresponding environment modes. They interact through beam splitters 
with transmittivity rj. If the modes are described by canonical operators qj,Pj, satisfying canonical commutation 
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FIG. 1: The model for a lossy bosonic channel C. Each input mode (left-right line), representing one use of the channel, 
interacts with the corresponding environment mode (top-bottom line) through a beam-splitter. To introduce memory effects, 
environment modes are initially considered in a correlated state. 
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For a memoryless channel, the environment modes are initially in an uncorrelated state, the vacuum state in the 
simplest case. More generally, we can consider a Gaussian state characterized by covariance matrix V env and zero 
displacement vector. In terms of Wigner function we have 
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where x env := (qf nv , . . . , q™ v , pl nv , . . . ,p n nv ) with denoting the transpose and | • | the determinant (when applied 
to matrices). Throughout the paper the normalization of a n-mode Wigner function reads J W(x) dx = (27r) n . 

Let us now consider the classical use of the bosonic channel where input quantum states carry the values of a 
random classical variable. Then, the mapping p \— > C[p] can be seen as a mapping of phase space points, so that for 
each mode we have two real random values (or one complex) carried. Therefore, we label an input state over n modes 
(channel's uses) as with a := (3?(ai), . . . , $t(a n ), 9(ai), . . . , 9(a„)) T and ay € C. 

Restricting our attention to input Gaussian states [g], we consider the classical variable a encoded via random 
displacements of a suitable (Gaussian) seed state <Ji n , that is 



P?n = [^Moj)] (Tin [®?=il>i(ai)] + 



where Dj denotes the displacement operator on the j-th mode and a is chosen according to a Gaussian distribution 
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: exp [-a T V d a] , 



(2) 



having classical covariance matrix V c i/2. 

In the memoryless case, the classical capacity is achieved by using a in = |0)(0|®" and V c i — diag(iV), where N 
represents the average number of input photons per mode Q . In such a case entangled inputs turn out to not be 
useful. 

Quite generally we shall consider a{ n characterized by a covariance matrix Vi n and V c i not necessarily be diagonal. 
Then, the state pf^ will be characterized by the Wigner function 



: exp 



■i (x in - \/2a) T Vr} (x in - V2a 



where x in := (q\ n , q™ ,p\ n 
by the Wigner function 



l ) T . In turn, the average input state, ~p in = J da P(a)p a i , will be characterized 



W(x in ) = / da P(a)W i ^(x in ), 
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having covariance matrix 

V in = V in + V cl . 

The output state p^ ut = C [pj%] is still Gaussian and it has covariance matrix 

V out = rjV in + (1 - rj)V en v, (3) 
as consequence of transformations ([T]). In fact, the Wigner function of the input and environment is given by 

1 



with block covariance matrix 



■. exp 



~ (*tot - V2a tot ^j V to l (x t ot - VZatot 



(4) 



Vtot := 



V in 

V c ,r,, 



and vectors x.J ot := (x,^, xj nv ), otJ ot := (ot T , T ). By applying the beam splitter unitary transformation x tot — > -Bx t 
in Eq.Q with the 4n x An block matrix 



and then integrating over the environment variables x e „„ , we arrive at the output state Wigner function 

WZMout) = J- — = exp -~ (x out - ^pfr) a) V m \ (x out - a 
y\v ou t\ L z v y v 

Finally, the average output state ~p out = J da P(a)p a lt = C [pi n ], will be characterized by the Wigner function 



(5) 



W ou t(x out ) = / daP(a)W£ tt {-x. out ), 



having covariance matrix 



Vout = Vout +r]V c i- 



(6) 
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III. THE MEMORY MODEL 



To introduce the memory effect among different channel uses, we consider the environment modes initially in a 
correlated state, actually a multimode squeezed vacuum. We are going to use multimode squeezed states depending 
on a single parameter as introduced in Ref . , so that we choose the covariance matrix V env as 



V, 



Here, the n x n matrix Q is taken to be 



exp(sfi) 
exp(— sfi) 



(8) 



/0 1 .. 
1 1 



n = 



1 o 1 



\o o 



o \ 






1 J 



(9) 



and the parameter set represents the memory strength (for s = we recover the memoryless case). Notice that 
the memory effect is symmetric among all modes. Moreover, it decays over the number of uses making the model 
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quite realistic. Actually, if we consider the set of V env matrix elements belonging to a fixed row as a discrete function, 
this is well fitted by a Gaussian having width 2^/[s[. Such exponential decay of correlations makes the noise in the 
channel non-Markovian and the channel not forgetful [Tl| . 

Since we would explore the usefulness of entangled inputs in the presence of memory, we then consider <ii n having 
covariance matrix like that of Eq.©, that is 



exp(r^) 
exp(— rfi) 



(10) 



where r £ 1 is the entanglement parameter (for r = we recover separable vacuum inputs). 

To respect the input energy constraint, the added number of photons per mode by entanglement, 



^Tr(F„ 
2n 
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must necessarily be a fraction of N, say 9 n N with < 6 n < 1. This limits the range of r to a finite interval of . 
particular, defining 

Tr(F in ) - n 
2nN 

the roots of the equation 6 n = 1 determine the minimum (negative) and maximum (positive) allowed values of r. 

Besides entanglement in input, we could introduce classical correlations as well through V c i [hence P{ot)\. Since 
entanglement adds a fraction 9 n N of photons per mode, the diagonal of Vd must contain a fraction (1 — 9 n )N of 
photons, while off diagonal could introduce correlations. 

We choose the covariance matrix V c i describing classical noise to be of the same form of V env and V^„ , that is 



where 



V P 1 



Y = - 
2 



= 2nN{\ - 9 n ) 
Tr(y) 

exp(yfi) 



Y, 



(11) 







exp(— yfl) 



with the parameter y EM. describing classical correlations. 

It is worth remarking that the choice of the form of covariance matrices V env , Vi n , and V c i allows us to diagonalise 
them in the same basis. The eigenvalues of this kind of matrices are shown in Appendix A. Moreover, in Appendix 
B we show that the amount of average photon number per mode remains finite in these matrices by increasing their 
dimension. 

IV. THE HOLEVO-x QUANTITY 

Recall that the von Neumann entropy of a Gaussian state p with covariance matrix V, can be simply written in 
terms of symplectic eigenvalues of V . This means that 



%)=X>(k 
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where 



and are the symplectic eigenvalues of V, i.e. ±«^ v ; are the eigenvalues of £ 1 V with £ the symplectic 2n x 2n 
block matrix 



g (x) := (x + 1) log 2 (x + 1) - x log 2 x, 
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The classical capacity of a memoryless quantum channel is determined by means of the output Holevo-x quantity 

Xn := S (p out ) - J d<xP( a )S (pZt) ■ (12) 

In case of memoryless quantum channels, this quantity represents the highest achievable rate for a fixed encoding 
scheme. Thus, its maximization overall possible encodings (i.e. possible input states and probability distributions- 
Gaussians for the conjecture of Ref.[l]) gives the classical capacity. 

Assuming that v plays the same role for memory channels as well and assuming the possibility to extend the 
conjecture of Ref.Q], we should maximize Eq. (fT2")) overall Gaussian states <7j„ and distributions P(a). As a first step, 
here we consider the maximization overall Gaussian inputs of the form of Eqs. pop and pip . 

Since the covariance matrices (J3|), ([7]) do not depend on a, we get \ simply depending on the two parameters r,y 
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the symplectic eigenvalues of V out and V ou t respectively. 



with i/. 

By applying the model of Section IHII we get the following result for the symplectic eigenvalues (see Appendix A) 
1 
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From a geometrical point of view, each symplectic eigenvalue of Eq. (|T4]) can be interpreted as the relation among 

triangle sides of length 77, 1 — 77 and 2v^f ^ with the angle between sides of length 7/ and 1 — 77 equal to 7r — itpifi . 
Analagously, each symplectic eigenvalue of Eq. (TT5|) can be interpreted as the relation among tetrahedron faces of area 



77, 1 — ?y, ?7AT„ and 2v\f ^ with dihedral angles iipifi , 
We can now define the quantity 
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and iif'sy among the first three faces. 



(18) 



C„ := max ~Xn(r,y). 

r,y£M 77 

which would be a marker of the behavior of the capacity for any number 77 of channel uses. Substituting Eqs.(Q2 
and (|15[) in the relation (|13[) and maximizing it over possible values of variables r, y we can evaluate the quantity Q18p . 
Such a quantity is plotted in FigJ^l-left versus 77 for different values of s. On the one hand, for 77 = 1 (or equivalcntly 
for s — 0) we recover the value g(r]N) of the memoryless channel capacity [9J. On the other hand, it is clear the 
possibility to overcome the memoryless channel capacity (by using entangled inputs) when s ^ 0. In such cases we 
may notice different transient behavior by increasing 77. 

By taking the limit 77 —> 00 of C n we find the asymptotic behavior. Actually, this limit results the limit of Riemann 
sums thus leading to an integral (as described in Appendix B). Due to the regularity of the function in Eq. (|13D we 
can write 
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where i/ y01 ") and \A V -* are given by Eqs. (TTJJ and fl~5|) with K n replaced by K (see Appendix B) and ip^ replaced 



by ipim ■■=2(1- m) cos£. 

In Figl^l-right the quantity lim n _ 



[max ye R -^Xn( r ,y)] is plotted versus r. The maxima correspond to values of 



the quantity (|19p . For s ^ 0, it clearly shows the possibility to overcome the memoryless channel capacity by using 
entangled inputs (r ^ 0). 
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FIG. 2: On the left, the quantity C n is plotted versus n. For each n, bars from left to right refer to s = 0, 0.8, 1.6, 2.5 
respectively. The other parameters are TV = 8 and r) = 0.7. On the right the quantity lim n ^oo [max^gu ^Xn( r i2/)] i s plotted 
versus r. Solid, dotted, dashed and dashed-dotted curves refer to s = 0, 0.8, 1.6, 2.5 respectively. The curves attain maxima 
at 4.05, 4.07, 4.12, 4 respectively. The other parameters are N = 8 and r\ = 0.7. 

V. CONVENTIONAL DECODINGS 

We are now going to consider information transmission rates by conventional decoding procedures like heterodyne 
and homodyne measurements @, fl2| . 

A. Decoding by joint quadratures measurements 

Let us consider decoding at the output by joint measurement of canonical variables on each mode (heterodyne 
measurement). This is described by the probability operator measure (POM) <8" =1 (I0)(0l/ 7r ) with \Q) coherent 
state of the j-th mode having complex amplitude Q. 

In analogy with the vector a we are going to introduce the vector £ = (3?(Ci), ■ ■ ■ , 3?(Cn), ^(Ci): • • • j ^(Cn)) T - The 
probability for coherent amplitudes C is related to the probability for the corresponding quadratures x out , which is 
given by the Wigner function |5]). Thus, the conditional probability of getting £ at the output given the encoded a 
at input results 



p(CH 



: exp 



(20) 



where 



V c \ a := V out + -, 

is the relation for quadratures' covariance matrices. Here, the term 1/2 represents the quadrature vacuum noise added 
by heterodyne measurement. Then, using Eqs.® and we get the output probability 
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The Shannon differential entropy for a stochastic variable <I> taking real values <f> distributed with probability P(<f>) 
reads 
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FIG. 3: On the left, the quantity F n (heterodyne case) is plotted versus n. For each n, bars from left to right refer to s — 0, 0.8, 
1.6, 2.5 respectively. The other parameters are N = 8 and r\ = 0.7. On the right the quantity linin^oc, [max ag K -J(Z : A)(r, y)] 
is plotted versus r. Solid, dotted, dashed and dashed-dotted curves refer to a = 0, 0.8, 1.6, 2.5 respectively. The curves attain 
maxima at 2.72, 2.64, 2.43, 2.24 respectively. The other parameters are N = 8 and rj = 0.7. 



Applying this definition to Eqs.([2]), (|2"0)) and (j2"Tj) we arrive at the mutual information [l£ 

I(Z : A) = H(Z) - H(Z\A) 

By using the model of Section [Till we explicitly get 
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I(Z:A) = -J2 Poga (F+, k ) + log 2 (F_, fe ) 
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The quantity F n := max r ye ^I(Z : A)/n is plotted in FigJHJ-left versus n for different values of s. On the one hand, 
for n — 1 (or equivalently for s = 0) we recover the value log 2 (l + f]N) of the memory less channel rate @. On the 
other hand, it is clear that in this case the memory effects worsen the rate with respect to the memoryless case. 

By taking the limit n — > oo of F n we find the asymptotic behavior of the rate. Actually, this limit results the limit 
of Riemann sums thus leading to an integral (as described in Appendix B). Due to the regularity of the function in 
Eq. ([24| we can write 
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In FigJHJ-right the quantity limn^oo [maxj, e K : A)(r, j/)] is plotted versus r for different values of s. The attained 

maxima correspond to values of F . We see that the maximum for s = implies r = (and also y — 0) [this is a global 
maximum for the mutual information]. Nevertheless, for s ^ 0, the maxima are obtained for r ^ 0. Hence, entangled 
inputs give a higher rate with respect to separable ones in presence of memory, although such a rate is always smaller 
than the one for memoryless case. 



B. Decoding by single quadratures measurements 



Let us consider decoding at the output by measurement of a single canonical variable on each mode. This is 
described, e.g. by the POM ®" = i (|3?(C?)) (3^(0)1)- I n the following it will be useful to write a = (qh, aj), and 
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FIG. 4: On the left, the quantity F n (homodyne case) is plotted versus n. For each n, bars from left to right re- 
fer to s = 0, 0.8, 1.6, 2.5 respectively. The other parameters are N — 8 and r\ — 0.7. On the right the quantity 
lim n _, 00 [maxy gK i/(5RZ : KA)(r, y)] is plotted versus r. Solid, dotted, dashed and dashed-dotted curves refer to s = 0, 
0.8, 1.6, 2.5 respectively. The curves attain maxima at 2.274, 2.27, 2.25, 2.21 respectively. The other parameters are N — 8 
and r) = 0.7. 



C = (CriCi)- Then, we have to consider information only encoded in a R . Hence, while keeping the usual energy 
constrain, we will consider the matrix Y in Eq. (|lip now written as 



Y 



e yn 




Interpreting the Wigner function {5| as the probability for £ conditioned to a, and taking into account that the 
stochastic variables an, £ R are independent from ctj, £ /; we get, by integrating Eq.((5]) over £ 7 , 



P(C R \a R ) 



t/2 



v {1 } ] 
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: exp 



f(v^) (Cr-Vv 

olr) 



(26) 



Analogously we can consider the Wigner function © as the probability for £ and by integrating it over £ 7 we get 



P(Cr) = 



r n/2 
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<»R y V out J V.J 
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Finally using the probabilities ([26]) and (|2T[) in Eq.([22|). we arrive at 

J(3tZ:3L4) = H{ftZ) - H (3tZ|3iA) 



By using the model of Section [TTT| we explicitly get 



IQftZ : SL4) = ^J2 lo §2 1 + 2K r. 
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The quantity F n :— maxj- ygg I(JftZ : 3?j4)/n is plotted in FigHl-left versus n for different values of s. On the one 
hand, for n = 1 (or equivalently for s = 0) we recover the value (1/2) log 2 (l + 4^^) of the memoryless channel rate Q. 
On the other hand, it is clear that in this case the memory effects worsen the rate with respect to the memoryless 
case. 

By taking the limit n — > 00 of F n we find the asymptotic behavior of the rate. Actually, this limit results the limit 
of Riemann sums thus leading to an integral (as described in Appendix B). Due to the regularity of the function in 
Eq. ([HI we can write 



F := lim F n — max — 



rj e 



2y cos £ 



cos « + (1 - rj) e 
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FIG. 5: On the left, the optimal r value for the quantities of Eqs. (|19p (solid line), (|25[1 (dotted line) and (|29p (dashed-dotted 
line) is shown versus s. The value of the other parameters are N — 8, rj — 0.7. On the right, the quantity C of Eq. (|19[) is 
plotted versus s for values of rj going from 0.1 (bottom curve) to 1 (top curve) with step 0.1. The value of the other parameter 
is N = 8. 

In FigdJ-right the quantity lim^oo [max^gR i/(3?Z : SRA)(r, y)] is plotted versus r for different values of s. The 
attained maxima correspond to values of F. We see that the maximum for s = implies r = (and also y = 0) 
[this is a global maximum for the mutual information as in heterodyne case]. Nevertheless, for s ^ 0, the maxima are 
obtained for r ^ 0. Hence, entangled inputs give a higher rate with respect to separable ones in presence of memory, 
although such rate is always smaller than the rate of heterodyne detection (for the same values of s) . 

VI. CONCLUSION 

We have presented a realistic model to describe memory effects in a lossy bosonic channel. These effects have 
nontrivial long range correlations (non-Markovian) . Notwithstanding, the model has allowed us to characterize the 
channel over an arbitrary number of uses for classical information transmission. Therefore, the asymptotic behavior 
of a Gaussian memory channel has been studied for the first time. 

We have shown the usefulness of entangled inputs in the presence of memory. In particular, we have shown that 
they can enhance the classical capacity with respect to the memoryless case. This has been done by maximizing 
the Holevo-x quantity over a class of Gaussian states. We have also shown the utility of entangled inputs for rates 
achievable by conventional decoding procedures (heterodyne and homodyne measurements), although these do not 
exceed the classical capacity of the memoryless channel. 

The usefulness of entangled inputs is summarized in FigHj-left where we can see that the optimal entanglement 
degree smoothly varies with the degree of memory. The absence of any kink excludes phenomena similar to phase 
transitions in contrast to what happen for qubit memory channels [l3| . Still, in FigfSJ-right we can see the enhancement 
of the classical capacity in the presence of memory for channels characterized by different values of transmittivity 
rj. Such enhancement grows thinner and thinner as r\ approaches one. Actually, the maximum in each curve shifts 
towards s — as rj approaches one. For rj — 1 the quantity C does not longer depend on s because the environment 
no longer affects inputs. We also notice that for increasing values of s, curves corresponding to different values of r\ 
(except -q — 1) flow together towards a single value. This is presumably due to the fact that for large values of s the 
class of Gaussian states used for the maximization is too small to get the optimal solution (which seems in agreement 
with the saturation of the optimal r value). 

In conclusion, we believe that this study paves the way for a deeper characterization of lossy bosonic memory 
channels. The next step would be the full maximization of the \ quantity overall Gaussian inputs. It remains an open 
question what could be the optimal decoding procedure able to attain a capacity higher than the memoryless one. 
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APPENDIX A 



Let us consider the n x n matrix O of Eq.©. To obtain its eigenvalues we first notice that Q, = 21 — T, where T is 
the following n x n matrix 



/ 2 -1 

-1 2 -1 ... 
0-12-10 



T := 








... 0-12-10 

0-12-1 

\0 -1 2 ) 



Being a finite-difference counterpart of the differential operator —-j-21 the matrix T plays an important role in 
computational mathematics and its properties are well known. For instance, it can be shown that it has the following 
eigenvalues : 



\P =2(1- 



n+1 



with corresponding normalized eigenvectors with components 



(T) 

v j,k 



2 . jkn 

■sin r, k = 1, ...,n. 



n+1 n + 1' 



(T) 

Thus, the matrix f2 will have the same eigenvectors v ) and its eigenvalues turn out to be 



Af ) =2-Af ) =2cos^. 
3 J n+1 



(30) 



For any matrix like 



A =2 



1 fe^ n 



e 



(31) 



we get the eigenvalues 



(A) = 1 e± 2 7 cos(^) 

±,k 2 



(32) 



as direct consequence of Eq. ([3D]) . 

Since V env , Vi n and V c i have the same block structure in terms of Q, they can be diagonalised in the same basis. 
Hence, by taking into account Eqs.©, ([7]), we straightforwardly get the eigenvalues of V out and V out 



*±5r = i(r ? e ±2rcos (^) + (l-7 7 ) e ±2scos (^) 



with K n given by Eq. (fT( 

Finally, by definition of symplectic eigenvalues we get 



(v out ) 

k 

(Vout) 



'k 



,(V out ) x (v out ) 

A +,k A -,k J 

A +,k A -,k ) 



leading to Eqs.(fl4]) and (fT5 
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APPENDIX B 



We show that dealing with matrices of the form of Eq. (|3ip , the average number of photons per mode remains finite 
even in the limit n — > oo. To this end we consider 

Tr(A) , Tr(A) 1/ Trfe^ 2 ) , Tr( e -^)\ . , 

lim — ±—L = lim —. v + = - lim v - ; + lim — > '- . (33) 

n— too 2n n^oo 2(71 +1) 4 yn-^oo 71+1 n^oo 71+1 / 

By taking into account Eq. (|3li|) we can rewrite the first term at right hand side of Eq. (f3"3")) as 

hm — = hm ^ fc - i 



n + 1 n— >oo 71+1 

This relation is the limit of Riemann sum becoming an integral for the function f(x) — e 2 7cos7rx This l eac \s to a 
modified Bessel function of the first kind and zero-order [15[ 

lim — i '- = / e 2 ~* cos ™dx = - / e 2 T cos «^ = / (2 7 ). 

n^oo 71 + 1 J TT J 

Since the Bessel finction I is even, the result of Eq. ([33|) is Iq(2^)/2 whose asymptotic behavior is e 2 ~* / (iy/ifj) for 
large 7. The existence of a finite limit in Eq. ([33|) allows us to conclude that the considered matrices (like (f3"Tj) ) give 
rise to a physical model. 

Finally, by applying the above method to calculate the limit values of 9 n and K n given by relations (| 1 T|) and (1 16|) , 
we get 

9 := hm 9 n = W- 1 , 

^M^J^^-W (34) 



n— »oo 



io(2tf) 
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